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Abstract. We report on the analysis of new X-ray data obtained with XMM-Newton and Chandra from two ROSAF-discovered 
X-ray dim isolated neutron stars (XDINs). RX J0806.4— 4123 was observed with XMM-Newton in April 2003, 2.5 years after 
the first observation. Fhe EPIC-pn data confirm that this object is an X-ray pulsar with 1 1.371 s neutron star spin period. Fhe 
X-ray spectrum is consistent with absorbed black-body emission with a temperature kF = 96 eV and Nh = 4-10 19 cm" 2 without 
significant changes between the two observations. Four XMM-Newton observations of RX J0420.0— 5022 between December 
2002 and July 2003 did not confirm the 22.7 s pulsations originally indicated in ROSAF data, but clearly reveal a 3.453 s period. 
A fit to the X-ray spectrum using an absorbed black-body model yields kF = 45 eV, the lowest value found from the small group 
of XDINs and Nh = 1.0-10 20 cm -2 . Including a broad absorption line improves the quality of the spectral fits considerably for 
both objects and may indicate the presence of absorption features similar to those reported from RBS1223, RX J 1605. 3+3249 
and RX J0720.4— 3125. For both targets we derive accurate X-ray positions from the Chandra data and present an optical 
counterpart candidate for RX J0420.0— 5022 with B = 26.6±0.3 mag from VLF imaging. 

Key words. X-rays: stars - stars: neutron - stars: magnetic fields - stars: individual: RX J0420.0— 5022, RX J0806.4— 4123 



^ ■ 1 . Introduction observation as part of the telescope scientist guaranteed time 



revealed a possible candidate for the neutron star spin period 



b The two soft X - ni > sources KX. 10X06.4 4123 and of ] 1 .37 s,' formally detecte d at a 3.5a level in the F.PIC-pn 
, RX J0420.0-5022 belong to the group of X-ray dim iso- data (Haberl & Zavlin 2002 hereafter HZ02) 
lated neutron stars (XDINs) discovered in the ROSAF all-sky 
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i leviews seei ireveseiai^uuut imoKa xhe faintest object in x . rays the XMNs ifj 

RXJ0806.4-4123 was discovered by rx J04 20.0-5022, which was selected initially due to a mis- 

a dedicated search for XDINs apply- identification with a nearby galaxy . FoUow . up ROSAT PSPC 

ing spectral hardness ratio selection criteria to the ROSAF and Hpj observations and ESO-NFF images, which revealed 

survey source catalogue restricted to the galactic plane. Fwo nQ Qptical counterpart brighter than B = 25 . 25i mled out possi . 

short ROSAF pointed PSPC observations covered the source We kinds f known X-ray em itters other than an isolated neu- 

serendipitously at an off-axis angle of 31 and the obtained ^ gtar (Haberl gt aJ lgg ^ The x . ray spectrum . consistent 

spectra were consistent with a black-body model with tern- wkh a planckian shape as observed from the other XDINs - ob- 

perature kT = 78±7 eV. No optical counterpart was identified tained from ^ short pspc pointed observation allowed only 

with a limiting B magnitude of -24. A first XMM-Newton CQarse estimates f b i ack - b ody temperature and interstellar ab- 

~Z , rr 7Z ~ T IT , . , . sorption, but indicated that it is one of the XDINs with low- 
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* Based on observations with XMM-Newton, an ESA Science est temperature. Fhe PSPC data also showed evidence for pe- 

Mission with instruments and contributions directly funded by ESA rlodlc modulation of the X-ray flux with a period of 22.7 s 

Member states and the USA (NASA). Optical observations were per- suggesting that RX J0420.0-5022 is a X-ray pulsar similar to 

formed at the European Southern Observatory, ESO programme 66.D- RX J0720.4-3 125 faaberl et all 19971) . the only pulsar known 

0128(A). among the XDINs at that time. 
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Table 1. Chandra ACIS and XMM-Newton EPIC observations. 



Target 




Instrument 


Read-out 


Filter 


Sat. Revol./ 


Observation 




Exp. 








IVlOQc 




UDS.-1U 


Date 


Start 


E n( i /TTT^ 


|KSJ 


RXJ0806.4- 


-4123 


EPIC-MOS1/2 


FF, 2.6 s 


Thin 


168 


Nov. 8, 2000 


13:43 


18:47 


18.0 






EPIC-pn 


FF, 73 ms 


Thin 






14:24 


18:53 


15.6 
















13:35 


18:54 


19.1 






APTQ T 


EE Q O c 

rr, j.z s 




DVUZjy 


Feb. 21, 2002 


03:33 


UO.JJ 


1 1 Q 






EPIC -MOS 1/2 


FF, 2.6 s 


Thin 


618 


Anr ?4 7003 


14:18 


21:10 


24.6 






r,r lL.-pn 


FF, 73 ms 


i nin 






14:44 


o 1 - 1 n 

Zl . 1U 


zz. / 






RGS 1/2 


Spectro 








13:49 


21:11 


26.4 


RXJ0420.0- 


-5022 


ACIS-S 


FF, 3.2 s 


— 


500238 


NOV. 11, 2UU2 


21:38 


03:47 


19.6 






EPIC-MOS1/2 


FF, 2.6 s 


Thin 


560 


Dec. 30, 2002 


03:39 


09:43 


21.7 






EPIC-pn 


FF, 73 ms 


Thin 






04:01 


09:43 


20.0 






EPIC-MOS1/2 


FF, 2.6 s 


Thin 


561 


Dec. 31, 2002 


21:55 


03:59 


21.7 






EPIC-pn 


FF, 73 ms 


Thin 






22:17 


03:59 


20.0 






EPIC-MOS1/2 


FF, 2.6 s 


Thin 


570 


Ian. 19, 2003 


16:42 


22:55 


22.2 






EPIC-pn 


FF, 73 ms 


Thin 






17:05 


22:55 


20.5 






EPIC-MOS1/2 


FF, 2.6 s 


Thin 


664 


Jul. 25, 2003 


21:22 


03:26 


21.7 






EPIC-pn 


FF, 73 ms 


Thin 






21:45 


03:26 


20.0 



Here we present the results from the analysis of new 
XMM-Newton and Chandra data from RX J0806.4— 4123 and 
RXJ0420. 0—5022. We derive accurate X-ray positions utiliz- 
ing the imaging capabilities of Chandra and derive X-ray spec- 
tral parameters from the high statistical quality spectra col- 
lected by the EPIC-pn instrument on board XMM-Newton. We 
discuss the measured properties in comparison with the other 
known XDINs. 



2. Chandra and XMM-Newton observations 

Using XMM-Newton Jjansen et alJ l200ll) we observed 
RXJ0806.4-4123 and RX J0420.0-5022 during the A02 
guest observer programme. Here we report on the analysis of 
the data collected with the European Ph oton I maging Camera 
based on a pn f EPIC-pn. Istruder et alJhoOll) CCD detector 
which is mou nted behind one of the three X-ray telescopes 
jAschenbach et aljEoOol) . We include the first XMM-Newton 
observation of RX J0806.4-4123 (HZ02) from the guaran- 
teed time programme for a comparative study. Spectra ob- 
tained from RX J0806.4-4123 (RX J0420.0-5022 is not suffi- 
ciently bright) by the Reflection Grating Spectrometers (RGS, 
Iden Herder et alJ boOlh are also investigated. Details of the 
XMM-Newton observations are summarized in Table \l\ We 
do not use the data from t he EPIC-MOS (EPIC-MOS1 and 
-MOS2, Turn er et alJ 12001 ) CCDs because of currently un- 
calibrated changes in the low-energy spectral response. With 
Chandra we observed the two neutron stars in Cycle 3 of 
the guest observer programme using the ACIS CCD arrays 
(Table[3. We used the ACIS data to derive precise X-ray posi- 
tions. Due to the low statistics (see below) we did not accumu- 
late energy spectra. 



2. 1. X-ray positions 

The best X-ray positions were derived from the Chandra 
ACIS pre-processed level 2 event files using the 'celldetect' 
task of the Chandra analysis package CIAO 3.0. To optimize 
the signal to noise ratio we accumulated images in the en- 
ergy bands 0.3-1.0 keV for RX J0806.4— 4123 and 0.3-0.7 
keV for RX J0420.0-5022 (see below). RX J0806.4-4123 
was detected with 986 net source counts (0.3-1.0 keV, cor- 
rected for point spread function losses) and we obtained 
RA = 08 h 06 m 23!40 and Dec = -41°22'30'/9 (J2000.0). For 
RXJ0420.0-5022 270 net source counts (0.3-0.7 keV) were 
found yielding a position RA = 04 h 20 m 01!95 and Dec = - 
50°22'48'.'l (J2000.0). The 90% confidenc e systematic posi- 
tion uncertainty is 0.6" ( Aldcroft et al. 2000). We compared the 
X-ray positions of four RXJ0420. 0—5022 field sources (de- 
tected at relatively large off-axis angles of 4' - 5.5') with coun- 
terparts from the USNO-B catalogue and found agreement to 
better than 0.5" in each coordinate. We did not find any system- 
atic shift which would allow to reduce the systematic position 
error. The Chandra position of RX J0806.4— 4123 is consistent 
with that derived from the XMM-Newton observation in Nov. 
2000 (HZ02). 

2.2. Average EPIC-pn X-ray spectra 

To obtain energy spectra with highest statistical quality we 
used the EPIC-pn data and processed the photon event files 
using a the XMM-Newton analysis package SAS version 6.0. 
Single-pixel events from circular regions around the sources 
with 30" (RX J0420.0— 5022) and 35" (RX J0806.4-4123) ra- 
dius were extracted and background was taken from similar 
nearby source-free areas. The SAS task 'arfgen' was used to 
calculate effective areas and ready-made redistribution matri- 
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Fig. 1. Merged EPIC-pn spectra of RX J0806.4-4123 and RX J0420.0-5022 accumulated from all available observations. Left: 
The best fit black-body model with photo-electric absorption is shown as histogram. Right: Model fit including an absorption 
line with Gaussian shape. The lower two panels in each case show the fit residuals on the same scale for each target (plotted in 
units of a and as ratio data/model). 



Table 2. Black-body model fits to the EPIC-pn spectra obtained from individual observations. 



Object 




Obs. 


Count Rate 


kT 


N H 


Flux 0. 1-2.4 keV 


X ?/dof 








[cts s -1 ] 


[eV] 


[10 20 cm~ 2 ] 


10~ 13 ergcrrT 2 s _1 




RXJ0806.4- 


-4123 


168-pn 


1.780±0.012 


95.6±1.1 


0.40±0.11 


28.8±0.2 








618-pn 


1.769±0.010 


95.5±0.9 




28.7±0.2 


1.50/273 


RXJ0420.0- 


-5022 


560-pn 


0.212±0.004 


44.8±1.8 


1.01±0.37 


4.8±0.3 








561-pn 


0.222±0.004 






5.1±0.3 








570-pn 


0.217±0.004 






5.1±0.3 








664-pn 


0.212±0.004 






4.9±0.3 


1.25/157 



Count rates with 68% errors correspond to the 0.13-1.2 keV (RX J0806.4— 4123) and 0.13-0.9 keV (RX J0420.0— 5022) bands as used for the 
spectral analysis. Fluxes are corrected for point spread function losses. Temperatures and fluxes are not corrected for relativistic effects and are 
given as measured by a distant observer. Errors on the spectral model parameters are derived for a 90% confidence level. 
"=" denotes fit parameter is a single free parameter common to all spectra of one object. 
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ces for EPIC-pn full-frame (FF) mode were used for the spec- 
tral analysis. 

2.2.1 . Black-body model 

The EPIC-pn spectra obtained from the two XMM-Newton 
observations of RX J0806.4— 4123 were fit simultaneously 
first with an absorbed black-body model. The column den- 
sity was treated as a free parameter common to both spectra. 
Temperature and normalization were allowed to vary freely 
for both spectra independently to allow for possible small gain 
variations (see also lHaberl et al.l2004l) and to investigate long- 
term flux variations. Both spectra are fully compatible in shape 
and absolute flux and the derived spectral parameters, which 
are summarized in Table [2] agree within the statistical errors 
(90% confidence for one parameter of interest). In compari- 
son to the fit presented by HZ02 the up-dated calibration yields 
reduced residuals and a more reliable value for the absorption 
(the Nh value derived with the same calibration from the EPIC- 
pn spectrum of RX J1856.4— 3754 agrees with that obtained 
from the Chandra LETGS spectrum). Because no significant 
differences are seen in the two spectra, we merged the data and 
accumulated a combined spectrum which is presented in Fig.Q 
Again it was modeled using the absorbed black-body spectrum 
and the best fit parameters are listed in Table[5] 

The four EPIC-pn spectra of RX J0420.0— 5022 were an- 
alyzed in a similar way as those of RX J0806.4— 4123. Due 
to the lower flux of this target we allowed only the black- 
body normalization to vary between the spectra of the indi- 
vidual observations and fitted temperature and Nh as com- 
mon parameters. The results are presented in Table [2] As for 
RX J0806.4— 4123 the four spectra are consistent with each 
other allowing to accumulate a total spectrum from the merged 
data (Fig.Ql. The fit results are listed in Table|3] The tempera- 
ture of 45 eV is the lowest value derived from any of the known 
XDINs. 

2.2.2. Absorption features ? 

The quality of the black-body fits to the phase-averaged 
merged EPIC-pn spectra of RXJ0806. 4-4123 and 
RXJ0420. 0—5022 is in both cases not very good. Recently, 
broad X-ray absorption features fro m the three XDINs 
RBS1223 (= 1 RXS J 1 30848 .6+2 1 270 8, lHaberl et all l2003t> . 
RX J0720.4-3125 jHaberl et alJl20()4 and RX J 1605. 3+3249 
Jvan Kerkwiik et al . 2004) were discovered in XMM-Newton 
spectra. In the following we add a broad absorption line to the 
black-body model. Because the statistical quality in the present 
spectra is relatively low we fixed the a width of the line at 70 
eV, a value similar to those found from RX J0720.4— 3125 and 
RXJ1605. 3+3249. The best fit parameters are summarized in 
Table |3] (phases 'all') and the model fit is presented in Fig.Q 
Formally, the quality of the fit improves considerably in both 
cases when including the absorption line. At the current stage 
of the EPIC-pn spectral calibration it is not clear how much of 
this improvement is caused by residual systematic calibration 
effects. In particular, the residuals of the black-body fit to the 



spectrum of RX J0806.4— 4123 are very similar in shape and 
strength to those of RX J 1856.4-3754 ( lHaberl et alJ2004 . On 
the other hand, for RXJ0420. 0—5022 the deviations are larger 
and lead to a stronger absorption feature (EW = —43 eV). At 
the derived energy of 329 eV RX J 1856.4-3754 does not show 
negative residuals which would lead to an artificial absorption 
feature. The reduction in the (reduced) \ 2 is also larger in the 
case of RX J0420.0-5022 than in RX J 1856.4-3754. 



2.2.3. Other models 

Following lHaberl et alJ (2003) we also tried non-magnetic 
neutron star atmosphe re models (e.g. iGansicke et alJ 12002: 
IZavlin & Pavlovl l2002h . As in the case of RBS1223 iron 
and solar abundance atmospheres do not reproduce the 
smooth continuum observed from RX J0806.4— 4123 and 
RXJ0420.0-5022 at energies above 0.5 keV where these 
models predict Fe absorption features. The Hydrogen atmo- 
sphere model fits yield for both sources lower temperatures as 
the black-body fits (35 eV for RX J0806.4— 4123 and 12 eV 
for RXJ0420.0-5022) which leads to th e well known prob - 
lem of predicting too high optical fluxes (Pav lov et alj fl996). 
The quality of the H atmosphere model fits range between 
those obtained from fits using a black-body with and without 
line (reduced \ 2 °f 1-7 for 137 dof and 1.5 for 51 dof for 
RX J0806.4— 4123 and RX J0420.0-5022, respectively). 

Another model, first applied by iBurwitz et alJ J2003I) to 
the spectrum of RX J1856.4— 3754 and more recently by 
Ide Vries et ail fc004 to spectra of RX J0720.4— 3125, consists 
of a black-body component attenuated by photo-electric ab- 
sorption and an additional ad hoc term oc E' 3 . This model can 
successfully describe the long-term spectral variations discov- 
ered in XMM-Newton RGS spectra from RX J0720.4-3125 
with the value of (3 increasing with time (over nearly 2.5 years) 
from 1.4 to 2.7. The reported value for RX J1856.4-3754 
of 1.28±0.3 is at the lower end of values measured from 
RXJ0720.4-3125. A fit of this model to the EPIC-pn spec- 
trum of RX J0806.4— 4123 results in an increased tempera- 
ture of 130 eV, an increased Nh of 1.9-10 20 cm~ 2 and = 
— 1.76, very different to the values reported from the high res- 
olution spectra of RX J1856.4-3754 and RX J0720.4— 3125. 
Formally, this model (with one more free degree of freedom) 
allows a better representation of the continuum curvature (re- 
duced x 2 = 1 -50) but does not improve the residual pattern seen 
on smaller energy scales in the pure black-body fit (Fig.^upper 
left). Hence, the best fit is obtained by including both, power- 
law attenuation and a narrower absorption feature (reduced \ 2 
of 1.0 for 134 dof). However, given the systematic calibration 
uncertainties it is not clear if the broad-band black-body model 
modification is justified in the case of RX J0806.4— 4123. For 
RXJ0420. 0—5022 the extremely soft spectrum covers a very 
narrow energy band which does not allow to constrain kT and 
P. Formally, the best fit yields (5 = -9.0 and kT = 550 eV. We 
therefore do not follow this model further for this source. 
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Object 




Instr. 


Phases 


Model 


kT 


N H 


Elinc 


EW 


X?/dof 












\C V J 


n 20 rm _2 l 


\C V J 


\C V J 




RXJ0806.4- 


-4123 


EPIC-pn 


all 


black-body 


95.6±0.9 


0.41±0.09 




- 


1.97/137 








all 


bb + gauss 1 


91.7±1.4 


1.12±0.22 


460±15 


-33±6 


1.51/135 






RGS 1/2 


all 


black-body 


104±4 


0.41 2 


— 


— 


0.91/245 








all 


1 1 i 

bb + gauss 


92±4 


1.12 2 


413±19 


-56±13 


0.82/243 


RX 10420.0- 


5022 


EPIC-pn 


all 


black-body 


44.9±1.3 


1.02±0.26 






1.85/51 








all 


bb + gauss 1 


45.0±2.6 


2.02±0.71 


329±24 


-43±3 


0.92/49 








soft 


black-body 


44.0±1.6 


1.02±0.32 














hard 


black-body 


45.7±1.6 








1.84/76 








soft 


bb + gauss 1 


44.5±1.3 


2.09±0.42 


329 2 


-45±10 










hard 


bb + gauss 1 


45.6±1.2 




329 2 


-46±8 


1.17/74 



1 Black-body model with Gaussian absorption line with fixed width of <r=70 eV. 2 Parameter is fixed in the fit. 



2.2.4. RGS spectra 

RGS spectra of RX J0806.4— 4123 were produced using 
'rgsproc' of SAS 6.0. The four spectra (RGS1 and RGS2 
from the two observations) were simultaneously fit (only al- 
lowing for individual normalization factors) using the black- 
body model with and without including a broad absorption line. 
Because of the lower sensitivity to photo-electric absorption 
in the RGS spectra Nh was fixed to the values obtained from 
the EPIC-pn spectra. The resulting best-fit parameters are pre- 
sented in Tableland the, for clearer representation combined, 
RGS spectrum is shown in Fig. 13 

Although the black-body model yields an acceptable fit to 
the RGS spectra, deviations are visible at low energies. The 
residuals are improved when the absorption line is added to the 
model. Different values for the line energy and the line strength 
are derived from RGS and EPIC-pn spectra, which most likely 
can be attributed to calibration problems. 

2.3. X-ray pulsations 

Due to the high time resolution of the EPIC-pn read-out mode 
only the data from this instrument were used for our tempo- 
ral analysis. We extracted pattern 0-12 events from circular 
regions around the sources with 30" (RXJ0420. 0—5022) and 
45" (RXJ0806.4-4123) radius. The event arrival times were 
checked for time jumps which sometimes occur in the pn data 
and corrected to solar system barycenter using the 'barycen' 
task. For p eriodicity searches w e followed HZ02 using the 
Z 2 method llBuccherietalJl983h and the Odds-ratio approach 
based on the Bayesian formalism ( Gregorv & Loredo 1996) for 
a precise determi nation of the spin period and its error (see 
IZavlin et all2004i for details). 

For both targets neutron star spin periods were reported 
in the past with low significance: 11.371 s at a 3.5er level 
for RXJ0806.4-4123 from EPIC-pn data (HZ02) and evi- 
dence for 22.7 s in ROSAT HRI data of RX J0420.0-5022 
llHaberletalJll999h . A search for pulsations around 11.371 
s in the data of the second XMM-Newton observation of 
RX J0806.4— 4123 confirms the pulse period and the Z\ peri- 



odogram and the pulse profile are presented in Fig. [3] We sum- 
marize in Table|4]the pulse periods with errors derived from the 
Odds-ratio method and the Z 2 values for RX J0806.4-4123 
with the values for the first observation taken from HZ02 for 
comparison. 

None of the four observations of RX J0420.0-5022 re- 
vealed a significant peak in the periodograms around the 
suggested ROSAT value of 22.7 s. However, in all four 
observations a period of 3.453 s was detected. As for 
RX J0806.4-4123, we performed Z 2 tests with contributions 
of higher harmonics. The results are listed in Table |4] with 
the best pulse periods and errors. The pulse periods are con- 
sistent within the la errors. The pulse profiles, folded at the 
best period and normalized to the mean count rates are shown 
in Fig. [4] Folding the light curves at the same average period 
yields no significantly different profiles. Although for all four 
observations the exposure, background level and net source in- 
tensity were nearly identical, the maximum Z 2 value varied 
from one observation to the next. However, the standard de- 
viation for the Z\ value expected from a sine-like signal is 
2 \fzf + 1, indicating that the variations are not significant. To 
verify this we performed an additional x 2 test following equa- 
tion (1) of lZavlin & Pavlovl fl998) by comparing the observa- 
tions pair by pair. Since the relative phasing is not known we 
stepped through a grid of phase shifts and used the minimum 
X 2 value to calculate the probability that the two light curves 
have different shapes. The x 2 values and inferred probabilities 
are summarized in Table|5] The highest probability (98.5% cor- 
responding to 2.4er), which was found from the comparison of 
the first and second observation, is still consistent with statisti- 
cal fluctuations. We conclude that there is no significant change 
in the pulse profile of RX J0420. 0—5022 over 0.5 years covered 
by the observations. 

The time baseline between the observations is not suf- 
ficiently long to detect significant period changes for both 
objects. Fitting a linear period change with time we de- 

s 

and 



rive spin period change rates 



and (0.5±8.6)-10- 



s s 



P of (-4.6±6.4)-10~ 12 
for RXJ0806.4-4123 



RXJ0420. 0—5022, respectively. The error on the spin period 
of RXJ0420. 0—5022 accumulates to nearly one period over 
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Fig. 2. RGS spectrum obtained by combining RGS1 and RGS2 spectra from the two observations of RX J0806.4— 4123. As in 
Fig.^the best fit black-body model with (right) and without (left) absorption line is plotted. 
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Fig. 3. Left: Periodogram derived from a Z\ test using the 0.12-1.2 keV EPIC-pn data of RX J0806.4-4123 obtained in satellite 
revolution 618 (April 2003). Right: Light curve folded at the best period, normalized to the average count rate. 



the 2 day gap between the first two observations. This is too 
large to be able to phase-connect the observations as it was 
originally one aim of the proposal where the gaps between the 
observations were optimized for a spin period of 22.7 s. 

2.4. Pulse phase resolved X-ray spectra 

For RXJ0420. 0—5022 no significant changes of pulse period 
and pulse shape were detected. Therefore, we merged the data 
of the four observations and performed a Z\ test to determine 
a medium period assuming P = 0. Light curves folded with 



the derived period of 3.453078 s were then produced for two 
different energy bands and a hardness ratio computed as ra- 
tio of the two light curves (hard band / soft band). These are 
shown in Fig. [5] (left). For the 0.12-0.7 keV light curve we de- 
rive (by fitting a sine wave to the pulse profile) a pulsed fraction 
of 12.7±1.8%. The pulsed fractions for the individual obser- 
vations are reported in the last column of Table [5] Variations 
in the hardness ratio (a x 2 test f° r a constant hardness re- 
sults in x 2 = 34.6 for 14 dof, corresponding to a 3.2cr result) 
are visible with the hardness ratio maximum delayed with re- 
spect to the intensity maximum, very similar to the behaviour 
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Fig. 4. EPIC-pn 0. 12-0.7 keV folded light curves of RX J0420. 0—5022 obtained from the four observations. The relative phases 
are arbitrary. 



Table 4. Results of the temporal analysis of EPIC-pn data. 



Table 5. Comparison of pulse profiles for RX J0420.0-5022. 



Obs. Period 
[s] 


zl 


zl 


7 2 


zl 


RXJ0806.4-4123 










168-pn 11.37139=b3.0-10 _4 


33.9 


42.9 


43.0 


43.2 


618-pn 11.37103±4.0-1CT 5 


43.6 


57.2 


60.4 


60.9 



RXJ0420.0-5022 

560- pn 3.453050±6.0-10" 5 30.0 36.9 48.5 51.7 

561- pn 3.453098±7.2-10" 5 27.2 30.1 32.1 33.9 
570-pn 3.453063±3.5-10~ 5 67.0 69.0 70.3 72.1 
664-pn 3.453074±8.4-10~ 5 36.4 36.5 39.0 41.5 

Errors on the period are given for a 68% confidence level. 



of RX J0720.4-3 125 ( Cropper et all200lllHaberl et alfcool) . 
The EPIC observations of RX J0420.0— 5022 cover only 0.5 
years and do not provide sufficient statistics to detect changes 
in the relative phasing between flux and hardness rat io, as re- 
ported for RX J0720.4-3 125 bv lde Vries et"ai1 ll2004 . 

To further investigate the spectral variations with pulse 
phase we accumulated spectra of RXJ0420. 0—5022 from 



Revolution 


560 


561 


570 


664 


PF [%] 


560 




21.97 


17.57 


10.99 


10.2zb3.5 


561 


98.5% 




3.68 


14.36 


10.0±3.5 


570 


93.7% 


3.9% 




8.61 


17.0±3.5 


664 


64.2% 


84.3% 


43.1% 




13.8±3.5 



The upper right part of the table gives the \ 2 values derived from the 
comparison of the two pulse profiles for 10 degrees of freedom, dof 
(12 phase bins with two parameters for the normalization). The corre- 
sponding probability that the two light curves of a pair have different 
shapes is given in the lower left part. 



phases around hardness maximum and minimum (phases 0.6- 
1.1 and 0.1-0.6 as used in Fig.|5}. The two spectra were again 
fitted with absorbed black-body model with and without ab- 
sorption line. The results are summarized in Table[5] As in the 
case of the phase-averaged spectrum adding the absorption line 
strongly improves the fit. There is an indication for an increased 
temperature during the phases of higher hardness ratio, but the 
errors are large and more statistics is required to prove the sig- 
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Fig. 5. Folded light curves in the total energy band covered by 
ratio, the ratio of the count rates in the hard to the soft band. 

nificance of this result. In Fig.[6]the spectra are presented with 
the best fit model including an absorption line. 

In the corresponding analysis of the merged data of 
RXJ0806.4-4123 the light curves, folded on the de- 
rived period of 11.371076 s are shown in Fig. [5] (right). 
RX J0806.4-4123 exhibits the smallest pulsed fraction of the 
known XDIN pulsars of about 6 % which renders the detec- 
tion of significant hardness ratio variations difficult. The fit of 
constant hardness ratio results in x 2 — 14.5 for 11 dof which 
corresponds to a probability of 79.3% (~ 1.3cr). However, it 
should be noted that by assuming a constant pulse period one 
may smear out the pulse profiles somewhat if the pulsar ac- 
tually exhibits a small period change over the 2.5 years from 
one observation to the next. Because of the marginal variations 
in hardness ratio we do not perform a phase resolved spectral 
analysis for RX J0806.4-4123. 

3. Optical observations 

Deep optical CCD images of the field containing 
RXJ0420. 0—5022 were obtained under photometric con- 
ditions on November 22, 2000 with FORS 1 at UT1 of the ESO 
VLT. Three CCD images with individual exposure times of 
1200 s were taken through a Bessel B filter and two images 
with 600 s exposure time each through a Bessel V filter. Two 
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the spectrum and in two sub-bands together with the hardness 
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Fig. 6. Phase resolved EPIC-pn spectra of RX J0420.0-5022 
fitted with a black-body plus Gaussian absorption line model. 

bright stars in the field have been masked by movable slitlets 
to avoid saturation of the detector and charge spilling. The 
mean seeing was 0.9" FWHM. 

The images were reduced and averaged using bias images 
and sky flats obtained during the same observing run and with 
the same instrument and detector setups. Cosmic ray events 
were removed by applying the FILTER/COSMIC command 
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Table 6. Optical photometry of candidate stars in the vicinity 
ofRXJ0420.0-5022. 



Object 


B [mag] 


V [mag] 


A 


24.40±0.04 


23.46±0.03 


D 


24.88±0.06 


24.38±0.07 


E 


25.33±0.10 


24.74±0.10 


H 


26. 11 ±0.20 


24.95±0.13 


I 


26.57±0.30 


>25.5 



of the ESO MIDAS software package to the averaged B- and 
V-filter images. An astrometrical calibration was derived from 
the positions of 21 field stars from the USNO-B catalog in the 
vicinity of RX J0420.0-5022. 

Within the error circle of the Chandra X-ray position of 
RX J0420. 0—5022 a single object is marginally detected on the 
B-filter image (Fig. |7J with B = 26.57±0.30 mag. No object 
was found at this position on the V-filter image down to a de- 
tection limit V > 2 5.5. Fluxes for objects A, D, E, and H from 
lHaberlet alj fl999). the new candidate star I, and the local sky 
background were determined using aperture photometry. For 
the photometric calibration nightly zero points and extinction 
coefficients were used from the FORS 1 quality control pages at 
ESO. The resulting B and V magnitudes are given in Table [6] 
The new B band photometry of stars A, D ; E, an d H is fully 
consistent with that given in Ha berl et al.l dl999l) . The detec- 
tion and photometry of star I is adversely affected by the lo- 
cally enhanced sky background due to the stray light halo of 
an 1 1 mag star at ~ 45" distance. The error uncertainty given 
for the B magnitude is the statistical error calculated from the 
photon statistics of the object and background fluxes. It does 
not include systematic uncertainties in the determination of the 
background flux which could be of similar order. 

4. Discussion 

XMM-Newton observations of the two isolated neutron 
stars RXJ0420.0-5022 and RX J0806.4— 4123 establish that 
both are X-ray pulsars with spin periods of 3.453 s and 
11.371 s, respectively. RX J0420.0-5022 exhibits a relatively 
strong sinusoidal intensity modulation o f ~12% similar to 
that observed fro m RX J0720.4-3125 llHaberl et alJ 1 19971 
IPaerels et all200ll) . Also, like for RX J0720.4-3 125, the hard- 
ness ratio shows variations with pulse phase with maximum 
(harde st spectrum) shifted with respect to the intensity max- 
imum dCropper et afll200ll ISaberl et alJEool Ide Vries et alJ 
2004). 

The light curve of RX J0806.4-4123 (Fig.QJ shows a pos- 
sible deviation from a purely sinusoidal shape, with a rel- 
atively smooth rise of the pulse followed by a steeper de- 
cline. The same behaviour may also be present (although to a 
lower level of confidence) in the light curve of RX J0420.0- 
5022 (Fig. |4}, and is different from what is exhibited by 
RX J0720.4— 3125. The pulse profile of the latter source ap- 
pears, in fact, q uite symmetric with respect to the pulse maxi- 
mum/minimum ( Haberl et alJl2004l) . An asymmetric (skewed) 




# * • i 

5 arcsec 

■ 1 * » w • ;5| 

Fig. 7. Averaged and smoothed B-filter image obtained with 
FORS 1 at the ESO VLT of the field of RX J0420.0-5022 with 
a total exposure of 1 h. The circle with 0.6" radius marks the 
Chandra position of the X-ray source and includes a possible 
optical counterpart with B = 26.57±0.30 mag. 



pulse profile hints towards a surface temperature distribution 
more complex than the (symmetric) one implied by a si mple 
dipolar magnetic field (Page ll995llPage & Sarmientoll996l) . If 
what we are observing is the thermal emission from the entire 
star surface, the simplest explanation would be to invoke the 
presence of higher order multipolar components in the mag- 
netic field of these objects. 

The pulse phase averaged X-ray spectra of both pulsars in- 
vestigated here are not very well fit by an absorbed Planckian 
model. The derived black-body temperature of 96 eV for 
RX J0806.4— 4123 is among the highest observed from this 
group of ROSAT discovered isolated neutron stars. The tem- 
perature kT of 45 eV inferred for RX J0420.0— 5022, which 
could only roughly be e stimated from the ROSAT PSPC spec- 
trum dHaberl e t aljl99S )). on th e other hand, is the lowest value 
observed to date Haberl ( 2004). Adding a Gaussian-shaped ab- 
sorption line to the continuum model improves the spectral fits 
for both objects, suggesting that there may be broad absorp- 
tion features presen t in the X-ray spectra similar to RBS1223 
dHaberl et alJ I2003I) . RX J160 5 3+3249 dvan Kerkwiik et alJ 
20041) and RX J0720.4-3125 dHaberl et all l2004 . If inteT 
preted as proton cyclotron absorption lines (as supported by 
the agreement between the magnetic field strengths indepen- 
dently inferred from spin-down measureme nt and cyclotron en- 
ergy in the case of RX J0720.4-3125, see lCropper et all2 004) 
the inferred energies of ^330 eV and ^460 eV would indi- 
cate magnetic field strengths of ~6.6-10 13 G and ~9.2-10 13 
G for RX J0420.0— 5022 and RX J0806.4-4123, respectively. 
However, because the residuals seen in the black-body fit to the 
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EPIC spectra of RX J0806.4— 4123 are similar in shape to those 
seen from RX J1856.4-3754 faaberl et al.l2004 and the RGS 
spectra yield somewhat different line parameters it is not clear 
to which extent they are influenced by systematic calibration 
effects. 

It is remarkable that probably up to five out of the six 
X-ray dim isolated neutron stars detected in the ROSAT all- 
sky survey, show broad absorption lines with central ener- 
gies in the range of ^200 eV (only an upper limit of ^300 
eV was found for RBS1223. iHaberl et alJl2003l) to 450 eV 
(RX J 1605. 3+3249. Ivan Kerkwiik et alJl2004 . An interpreta- 
tion in terms of fundamental lines of a proton cyclotron emis- 
sion implies surface magnetic fields in the relatively narrow 
range of B = 3-7 10 13 (l+z) G, with z ~ 0.3 the redshift pa- 
rameter. These estimates are comparable with the QED critical 
value Bqed — 4.4-10 13 G, that is expected to quench pulsar 
activity. The upper B value could simply reflect the lower sen- 
sitivity to line detection in the Wien part of the spectrum. Were 
these stars active pulsars, the spin periods of the order of 10 
s observed in some so urces imply t hat any radio pencil beam 
should be very narrow (lBi ggsll990l) and the chances then high 
that it does not sweep over the Earth. Magnetic field strengths 
below the critical value would still be compatible with the pres- 
ence of cyclotron lines if some of the detected features are 
higher order harmonics. Unfortunately, little is known about 
harmonics to fundamental line ratios. All opacities given so far 
for proton cyclotron lines have been computed neglecting the 
contribution of higher harmonics and the absence of such lines 
in model atmospheres does not mean that they should not be 
obser vable. In t he illustrati ve but s till unique case of IE 1207.4- 
5209 llBignami et al.l2003llPe Luca et alJ2004l) . the line equiv- 
alent widths of the various components (65 - 100 eV) do not 
vary much with central energy (690 to 2800 eV). Harmonic 
lines with similar equivalent widths could be missed in the very 
soft spectrum of RXJ0420. 0—5022, or maybe in the case of 
RX J 1605. 3+3249, with the highest line energy, but would have 
been detected in RBS1223 and RX J0720.4-3125. 

In the magnetic field and temperature conditions prevailing 
in XDINs, hydrogen atmospheres could be partly ionized as a 
resul t of the much increased i onization potential at large B (see, 
e.g., IZavlin & Pavlovll2002l and references therein). Bound- 
free and bound-bound transitions occur at energies which can 
be made comparable to those of the lines observed in XDINs. 
Since atomic line energies vary more slowly with magnetic 
field than cyclotron lines do, a larger range of surface mag- 
netic fields could in principle be allowed. The main constraint 
of this interpretation is that the proton cyclotron line should be 
located outside of the observable energy range, i.e., at E < 100 
eV, implying s urface magnetic fields below 210 13 G. Model 
atmospheres bv lHo et all J2003ll show that at B = 10 13 G, the 
s0-s2 transition at E ~ 180 eV and the H ionization edge at 
E ~ 250 eV could both contribute to some of the observed 
lines. At higher ene rgies, heavier elements could play a role 
jHailev & Moril2002l) . On the other hand, at such strong mag- 
netic fields (>3-10 13 G) and rather low surface temperatures 
(<10 6 K) the neutron star surface may be in a solid or con- 
densed state. At present, only a first step has been undertaken in 



studying radiative properties of such surfaces dZane et al.l2004t 
iTurolla et al.l2004l) . 

The precise X-ray position of RXJ0420. 0—5022 obtained 
with Chandra allowed us to identify a possible optical coun- 
terpart on a deep VLT B band image. The shorter V band ex- 
posure is not sensitive enough to provide colour information 
which could be used to support the identification. However, the 
extrapolation of the X-ray spectrum into the optical results in a 
factor of 12.3 (or smaller if our identification is wrong) lower 
flux than actually is observed in the B band. Similar optical 'ex- 
cess' in the spectral energy distribution is seen from all optical 
counterparts (and candidates) identified so far (factor 14 for 
RX J 1 605 .3+3249 ([Kaplan et al J l2003al) . <4.9 for RBS1223 
dKaplan et alJ 120021) . 2.4 with an additional c omponent in- 
creasin g with wavelength for RX J0720.4 — 3125 ijKa pian et all 
2003b) and -7 for RX J1856.4— 3754 JBurwitz etal.ll2003[ 
based on the spectral parameters derived from the LETGS X- 
ray spectrum). Understanding the overall spectral energy distri- 
bution of XDINs from radio to X-ray wavelength will be cru- 
cial to derive a comprehensive model of the emission properties 
if this group of neutron stars. 
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